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A Green's function approach is presented for the D-dimensional inverse square potential 
in quantum mechanics. This approach is implemented by the introduction of hyper- 
spherical coordinates and the use of a real-space regulator in the regularized version 
of the model. The application of Sturm-Liouville theory yields a closed expression for 
the radial energy Green's function. Finally, the equivalence with a recent path-integral 
treatment of the same problem is explicitly shown. 

1. Introduction 

The use of path integrals for the analysis of the bound-state and scattering sect( 
in quantum mechanics has significantly advanced in the past three decades. EH 
However, much work remains to be done, a situation that is particularly evident for 
the case of singular potentials. □ Among these, the contact interactions&ElEl and the 
inverse square potentialoB stand out; these problems have been the subject of ac- 
tive research in the past few years, mainly through the direct use of the Schrodinger 
equation combined with a pp ropriate regularization and renormalization analyses a 
la field theoryfflfflSffl Recently we have started a program centered on the 
use of path integrals for a systematic study of bound states and of the connection 
between field theory and quantum mechanics. EZnla We expect that our approach 
may shed light on the formidable problem posed by bound states in quantum field 
theory.E3 In addition to introducing a technique of infinite summations of perturba- 
tion theory, we have specifically found analytic solutions_for the two most outstand- 
ing singular interactions: Jihe inverse square potentialO and the two-dimensional 
delta-function interactionJla 

In this paper we now investigate the quantum-mechanical properties of the in- 
verse square potential from a complementary viewpoint, using an intrinsically non- 
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perturbative operator approach. This is accomplished by directly considering the 
problem associated with the energy Green's function G(E), defined as the Fourier 
transform of the path-integral propagator K(T) [multiplied by the Heaviside func- 
tion 0(T)]. The Gree n's fu n c tion is then computed by solving a Sturm-Liouville 
problem in real space£30'c3 Our final results are in complete agreement with 
the ones we obtained earlier by means of a path-integral perturbative approach.EJ 
More precisely, we explicitly show that both approaches reproduce identical energy 
Green's functions. 

The relevance of this computation is highlighted by its recent application to a 
simple realization of a auantum anomaly in Nature: the interaction of an electron 
with a polar molecule.Ej As shown in Ref. ^3|, the scale symmetry properties of the 
inverse square potential are inherited by the full-fledged dipole potential describing 
the physics of this system — electron binding occurs for a sufficiently strong coupling 
(dipole moment) and this amounts to quantum mechanical symmetry breaking. Our 
calculations in this paper provide additional support for this remarkable and unusual 
result. This mechanism is further analyzed in Sec. 4. 



2. General Framework 

The path-integral treatment for a D-dimensional nonrelativistic particle is based on 
the quantum-mechanical propagator 
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Equation (Q) has been explicitly formulated in Cartesian coordinates by means 
of a time lattice tj = t' + je [where e = it" — t')/N, with j = 0, • • • , N, while 
t"], such that Tj = r(tj), with the end points being ro = r' and 

N-l a(N) 



t = t' and t N - +" 



r^v = r". The corresponding discrete action in Eq. ([!]) is S^ N ^ — X^q 1 



with 
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M(rj + i~rj) 2 /2e — eV(rj,tj), for a particle of mass M subject to a potential 



V(r,t). The connection with the operator approach is established by means of 
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where T is the time-ordering operator and H is the Hamiltonian. For a time- 
independent Hamiltonian, the dependence of Eq. (|^) with respect to the times t' 
and t" is only through the time difference T = t" — t' . For this all- important case, 
the energy Green's function is then defined as a Fourier transform of the Green's 
function obtained from (0)-(H) by multiplication with the Heaviside function 9(T), 



i.e., 



i r°° 

G d (t",t';E) = 4z / dTe lET / h K D (r",ri;t",t>) 
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Equations (g) and (j3|) establish the operator form of the energy Green's function in 
its configuration-space representation, which explicitly becomes 



G D (r",v';E) 



E-H + ie 
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where ie = i0 + is a small imaginary part. The operator structure of Eq. (|4|) implies 
the differential equation 



| V 2 , + ™[E- V(r')]} Gn(r", r'; E) = 6{v" r') , 
in which we have conveniently defined the rescaled Green's function 



(5) 



G D (r",r';E) = —G D (r",r';E) 
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For the particular case of a central potential, Eq. (||) can be solved by separation 
of variables. This procedure can be systematically implemented by considering 
a complete set of angular functions, the D-dimensional hyperspherical harmonics 
Yi m (£l), for which I is the angular momentummrantum number and m = 0, . . . , di, 
with d { = (21 + D - 2) (I + D - 3)\/U(D - 2)!.EJ Then, the partial-wave expansion 

OC d[ 

G D (r",r';E) = (/y)-^^ £ £ Y lm (n")Y^(n')G l+1/ (r",r';E) (7) 

1=0 m=i 

implicitly defines the radial energy Green's function Gi+ U (r" ,r'; E) for each angu- 
lar momentum channel I and dimensionality D — 2(y + 1). Equation (0) r has been 
written to explicitly display the property of interdimensional dependences except 
for the scale prefactors (r'V) ^ D 1 ^ 2 , the only dependence of Gd{y" , r'; E) upon 
the dimensionality D is through the combination I + v, and this information is con- 
veyed by the function Gi +U {r" , r'; E). Likewise we define a rescaled radial Green's 
function Qi +y (r" , r'; £7) from the expansion of Gd{y" , r ']E), just as in (fjj) — or alter- 
natively, by directly enforcing the analogue of Eq. (^). Then, substitution of Eq. ([?]) 
in (|^) yields the radial differential equation 

+ f P - Wl - »±^} e , + ,(r»,r' ; E , . J( r» - /, , (8) 
which stands for a particular case of the Sturm-Liouville Green's function proble 

+ V + q(r')\g(r",r';E) =8(r"-r') , (9) 
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with p(r) = 1, 77 = 2ME/H 2 , and g(r) = -2MV c ff MZ^ 2 ; here KffO) includes, in 
addition to the potential V(r), the centrifugal barrier ,c3 displayed as the last term 
on the right-hand side of Eq. (N). 
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Equation (0) can be solved by applying the following standard technique. Let 
u l < u( r ) an d u iJ( r ) be the particular solutions to the corresponding homogeneous 
differential equation 



d 
dr 
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subject to the appropriate boundary conditions: Uj^(r) at the left boundary and 
u lv v) a ^ the right boundary. Then, the Green's function is given by 



Q l+v {r",r';E) 
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p(r>) W{u$, 



(11) 



where r< (r>) is the smaller (larger) of r' and r", while W |tt|^, ( r ) is the 

Wronskian of u^\r) and u 

Solving Eq. (^) requires specific knowledge of the potential and has to be dealt 
with on a case-by-case basis. We now turn our attention to the advertised problem: 
the inverse square potential. 

3. Green's Function for the Inverse Square Potential 

For the case of an inverse square potential, we conveniently write 
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As a particular case of the Sturm-Liouville problem (||), Eq. ( [i"3| ) can be solved 
by introducing the functions uf'J {r) and UiJ{r) that satisfy the corresponding 
homogeneous differential equation 



(l + v) -A - 1/4 



dr 2 



ui, v (r) = , 



(14) 



along with the boundary conditions outlined in Sec. 2. Equation (|l4| ) is immediately 
recognized to have solutions of the generic formil 



uiA r ) = VrZ st (y/rjr) 



(15) 



where Z Sl (z) represents an appropriate linear combination of Bessel functions of 
order 

Sl = A W -A , (16) 
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with 




inverse square potential V(r) = — h 2 X/2Mr 2 to a free particle is a shift in the 
angular momentum quantum number, 



This analysis shows that A, in Eq. ([17]) plays the role of a critical coupling, i.e., 
the nature of the solutions changes abruptly around the value A = Ar , for any 
state with angular momentum 1x3 Thus, x\*^ represents the threshold separating 
the two coupling-strength regimes: 

• subcritical or weak coupling, A < A; (i.e., real s/); 

• supercritical or strong coupling, A > A; (i.e., imaginary si). 

In the latter case the order can be rewritten as si = with a real parameter 



Classification of the nature of the solutions can be fully accomplished by also con- 
sidering the values of the energy E = H 2 rj/2M: 

(i) scattering states amount to r\ = k 2 > 0, for which the Bessel functions in 
Eq. (^5[ ) have real argument kr; 

(ii) bound states amount to rj = — k 2 < 0, for which the Bessel functions have 
imaginary argument kr = inr. 

As a result, the solutions ( |l5| ) fall into one of the following four families: 

(i) subcritical bound-state sector: 



I + v -> si = ^/(l + v) 2 - A . 



(18) 




(19) 



{l sl (Kr),K sl (nr)} ; 



(ii) subcritical scattering sector 



ui, v (r)/Vr 



(iii) supercritical bound-state sector: 



ui tV (r)/y/r 



{hsM T )-> K isM r )} ; 



(iv) supercritical scattering sector: 



{Hg\(kr),Hg\(kr)} . 
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Here the symbol {,} represents a linear combination and H^' 2 \z) stand for the 
Hankel functions of the first and second kinds, whereas I Sl (z) and K Sl (z) stand for 
the modified Bessel functions of the first and second kinds, respectively.t3 

The analysis of Ref . [l5] applies almost verbatim to the Green's function formula- 
tion, with the result that regularization is only needed in the supercritical coupling 
regime. However, in this paper — just as in Refs. [l2] and [l3| — we regularize the in- 
verse square potential by introducing a real-space regulator a, a procedure that 
is different from the dimensional regularization of Ref. [l5|. The introduction of a 
regulator is necessary for strong coupling because the inverse square potential fails 
to provide a discriminating boundary condition at the originJl3 Thus, this singu- 
lar behavior leads to the emergence of a scale — a behavior known as dimensional 
transmutation.ca In some sense, the key to a successful regularization procedure is 
the restoration of a sensible boundary condition at the origin. In particular, a real- 
space regulator allows one to restore the usual homogeneous boundary condition, 
but now displaced away from the origin to a radial position r = a. In short, the 
boundary conditions satisfied by the Green's function are: 

• boundary condition at infinity, 

i&'M^O; (20) 

• boundary condition at r = a, 

u ( <\a) = 0. (21) 

In what follows and as anticipated for the supercritical regime, &,= with 0/ real 
and given by Eq. (|l9|). From the behavior of Bessel functionscJ and the boundary 
conditions, one concludes that, for the bound-state sector of the theory, 

u}J> (r) cx ^ [K„ («a)7 ai («r) - I., (na)K Sl M] (22) 

[which satisfies condition (pi])] and 

v${r)ocViK ai (Kr) (23) 

[which satisfies condition (p0|)]. In addition, from the well-known properties of 
the Sturm-Liouville equation, the Wronskian of any two Bessel functions Z(z) is 
proportional to 1/z; for the problem at hand, it suffices to know (for example, from 
the small- argument behavior of the Bessel functions) that W {K s (z), I s (z)} = l/zx3 
Then, Eq. (||) yields straightforwardly the desired regularized radial energy Green's 
function 

/hi i . 2M i KieAnr^ 

G l+v (r",r';E X;a) = V^V 7 >; 

h z K i0l {Ka) 

x [Kie^Kajlie^Kr^ - I iei (Ka)K i&l («r<)] , (24) 
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with k = y/-2ME/h. 

Analysis of the scattering sector requires an analytic continuation to positive 
energies, with a replacement of the Bessel functions in Eq. ( pi[ ) as follows: I s {Kr) — > 
(-i) s J s (kr) and K s (nr) -> m' +1 HP(kr)/2; then, 



G l+v (r",r';E\X;a) 



2M TTl 



(ka)J ie , (Ar<) - J 40! {ka)H&> (fcr<) , (25) 



gg(^>) 



where fc = s/IME/h; from this, the S-matrix can be derived with the formulation 
of Ref. 

Equations (M) and ( p5| ) provide the solution to our problem because the energy 
Green's functions summarize all the physical information about bound states and 



scattering. Indeed, Eq. (24) is identical to the corresponding formula recently found 



in Ref. 17. This completes the proof of the equivalence between the perturbative 
method of Ref. [I?] and the nonperturbative approach of this paper. 

4. Renormalization, Anomalies, and Conclusions 

The analysis that follows from Eq. (g4|) is identical to that of Ref. [Tt], which we 
summarize here for completeness. The bound-state sector of the theory is deter- 
mined by the poles of Eq. (|24|), which provide the implicit solutions of the equation 
K SI (na) = 0. This equation provides no bound states in the subcritical regime, 
but the imaginary index si = iQi changes the nature of the Bessel function and 
leads to bound states in the supercritical regime. Then, as a is just a real-space 
regulator, attention needs to be focused only upon the limiting behavior for small a; 
consequently, the energy levels become E nr i = — (2e _7 /a) exp (— 2irn r /Qi) (with 
n = n r being the radial quantum number). This result, properly renormalized [with 
Qi = Qi(qX in the limit a — > 0], reproduces the known solution of the path-integral 
approachEII and of the Hamiltonian operator approach.li3iij'Ej 

A remarkable property of this solution is the emergence of an arbitrary dimen- 
sional scale, within a theory that is devoid of dimensional scales at the level of the 
Lagrangian — the phenomenon of dimensional transmutation.E3 

An even more illuminating approach to this scale problem is provided by a gen- 
eral analysis of symmetrvof theaction (but not necessarily the Lagrangian) under 



time reparametrizationsJ£3ti3 : E3a Invariance of the action leads to_the selection of 
the SO(2,ll conformal group, as for the inverse square potential,LJ the magnetic 
monopolej^j and the magnetic vortex.Lj The corresponding generators are: 



(i) the Hamiltonian H , associated with time translations t — > t 

(ii) the dilation generator 



a: 



D = tH 



-(r-p + p-r) , 



(26) 
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associated with a scale transformation t ~ * rt; r — > r 1 / 2 r; 
(iii) the special conformal generator 

if = i?t 2 - i(p • r + r • p)i+ ^mr 2 , (27) 

associated with the special conformal transformation with respect to time 
l/£-> 1/i + a. 

These three generators satisfy the commutator relations: [D, H] = —ihH, [D, K] — 
ihK, and [if, if] = 2ihD, which show that they form an SO(2,l) algebra. The 
same algebra is displayed by the two-dimensional delta-function intcraction.1 1 "Irl 
The anomalous version of this algebra as a consequence of the breakdown of the 
classical SO(2,l) symmetry is currently being investigated and will be reported 
elsewhere. 

This quantum symmetry breaking of SO(2,l) invariance is exhibited in Na- 
ture at the level of molecular physics, where it is realized in the point-dipole 
approximation.^ In particular, the interaction between an electron and a polar 
molecule can be reduced to a radial problem with an inverse square potential, for 
which the analysis of this paper applies. This information is then transferred to the 
angular part of the problem, which provides a critical dipole moment for electron 
capture and formation of anions — a result of unusual relevance and which is con- 
firmed by a large body of existing experimental and numerical evidence. EH This 
is indeed the most remarkable example of a quantum anomaly.EH 

In summary, we have completed a derivation of the energy Green's function for 
the regularized inverse square potential (using a regularization a la field theory) 
in the (strong) supercritical regime. This problem provides an example illustrating 
that the use of intrinsically nonperturbative techniques (operator Green's function 
approach) yields results identical to the corresponding ones derived from the infinite 
summation of perturbation theory to all orders. 
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